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The energy distributions of electrons of about 54, 75, and 92 

MeV were measured before and after passing through lead absorbers of 

2 

thicknesses ranging from 0.706 to 2.825 gm/cm , and through copper 

2 

and aluminum absorbers up to 5.726 gm/cm . The electrons were acceler- 
ated by the LINAC of the Naval Postgraduate School. A measurement to 
determine the optimum location of the absorber indicated that the 
separation between the scattering foil and absorber should not exceed 
6 cm, if geometric difficulties are to be avoided. The most probable 
energy losses agree with the theory of Blunck and Westphal for all 

three materials and all thicknesses. The half-widths of the distri- 

2 

butions agree with theory up to thicknesses of about 2 gm/cm . For 
greater thicknesses the experimental half-widths of copper absorbers 
appear to agree with theory; the experimental half-widths of aluminum 
are greater than theoretical predictions; and those of lead are smaller 
than theory predicts. 
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I. INTRODUCTION 



The energy distribution of an initially monoenerget ic electron 
beam which has passed through a layer of matter has been calculated by 
Blunck and Westphal [l]. This calculation evolved from the earlier 
theoretical work of Landau [2], EygerstS], Bethe and Heitler [4], and 
Blunck and Leisegang[5] . It essentially accounts for the energy loss 
distribution due both to ionization and excitation and to radiation. 

The theory is based on the assumption that the absorbing material is 
sufficiently thin so that the energy loss is small compared to the 
initial energy of the electrons. The development of the theory is 
contained in Section II. 

A number of experiments to measure the energy loss distribution 
of high energy electrons in thin layers of various materials have been 
performed, with incident energies in the range 10 MeV to 150 MeV. Of 
particular interest are the works of Breuer on aluminum at the 
Darmstadt linear accelerator [6]; Bumiller, Buskirk, Dyer, and Miller 
on aluminum at the LINAC of the Naval Postgraduate School [12,16], 
and Goodwin on copper, also at the Naval Postgraduate School [l3]. 

The experimental works listed above are generally consistent in their 
results . 

In each case agreement was found between experiment and theory for 

2 

thin absorbers (< 2.0 gm/cm ). However, for thicker layers of material 
the theory fails to predict the energy distribution experimentally 
realized. It is not clear how rapidly the theory breaks down as the 
energy loss becomes significantly large compared to incident energy. 



9 



In this thesis data are presented for the energy loss distribution in 

2 

lead for absorber thicknesses from 0.706 to 2.825 gm/cm . Additional 
data were obtained for thick layers of aluminum and copper to supple- 
ment the experimental results obtained by Bumiller et.al. [15] and 
Goodwin [13*]. The energy losses ranged from 17> to about 177> of the 
incident beam energy. The energy loss distribution can be character- 
ized by the most probable energy loss and the half-width, which is 
the energy width of the distribution curve measured at one-half the 
peak intensity. These values are used to compare experimental results 
with theoretical predictions, as well as to compare with other work. 
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II. THEORETICAL CONSIDERATIONS 



The Blunck and Westphal theory of the distribution for the energy 
loss suffered by an initially monoenerget ic beam of electrons in pass- 
ing through a layer of absorbing material assumes that the energy 
loss Q is small compared to the initial energy, E^. Let W(Q)dQ be the 
probability of energy loss between Q and Q + dQ and x equal that 
amount of Q lost by radiation. Then the ionization loss equals 
(Q - x). This leads to the equation 

c Q 

W(Q)dQ = ^ W ]; (Q - x)W s (x)dxdQ (1) 



where W_ and W are the energy loss distributions for ionization and 
I b 

radiation , respectively . 

For the distribution of energy loss due to ionization and excita- 
tion, the Landau equation [2], as modified by Blunck and Leisegang 

is used. The distribution is expressed as a function of the dimension 
less parameter X: 



W (Q)dQ = $ (X ) dX = 



r- c y 
\ n n 

njb 2 + Y* 



exp 



- V 2 ' 
2 2 
(b + Y„) J 



d\ 



( 2 ) 



where 

(3) 

The quantities introduced above are defined as follows: c , v , 

n n 

and X^ are constants used by Blunck and Leisegang to fit the Landau 
distribution to a sum of Gaussian functions. Their values are listed 
in Ref. [5} R is the thickness of the absorber measured in cm. 
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The quantity a is a function of the atomic number Z, the atomic 
weight A and the density p of the material, and (3 (= v/c) of the 
e lectrons : 



0.154 Z P MeV 

A (3 2 



(4) 



The quantity b is a correction to Landau theory given by Blunck 
and Leisegang [3"]: 



9 9 n V- 1 N 

, l 3.0 \ mm 

b = I r~ L~ 

m 



in 



2E. 



i (i - p ) 



(5) 



where the summation is made over the m ionization potentials of the 
atomic electrons of the material, while N is the number of electrons 



m 



with ionization potential 1^. Blunck and Westphal [l] and Breuer [6] 



used the approximation 



b = 2 x 10 



-5 QZ 



k„4/3 



(aR) 



( 6 ) 



This approximation is valid for aluminum and was used in making the 
theoretical calculations so as to permit comparison with previous ex- 
periments. However, the calculations for copper and lead were made 
using eq . (5) . 

Q is the average energy loss due to ionization (no radiation) for 
electrons of incident energy E^, and is given below. 

As can be seen from eq . (3), X depends directly on Q. Thus the 
width and most probable energy loss corresponding to the distribution 
$ (X ) depend on Q. Sternheimer [7,8,9] has calculated this quantity as 
follows : 



A t 




B + 0.43 + jjn E, - 0 2 
1 



c 



a s (x l 



m 



log 10 me-* 



MeV (7) 
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where t is the thickness in gm/cm and the constants As, B, C, X^, , 

and m g are parameters of the material and listed in Ref. 9. The most 
probable energy loss due to ionization and excitation (without radiation) 
is given by: 



A t 
Q . -2- 

» e 2 



B 



A t 
s 



+ 1.06 + 4n ^— —J - 



B - C • a lx 



Vi • log io —J 



MeV (8) 



For W , the distribution of energy loss due to radiation, Blunck and 



Westphal give [3,4]: 



H s (Q)dQ - 1 »r(.£) f 



(9) 



where 



a = 1.40 x 10 



-3 PZ 



, / 183\ 

2r Vl/3/ 



+ 1/9 



-1 



cm 



(10) 



and B is a normalizing factor = 



F(c^R + 1) 



The distribution of energy losses according to Blunck and Westphal 
is obtained by putting Eqs. (9) and (2) into eq . (1) and performing the 

integration. The results give an energy loss probability distribution 
for a single electron of incident energy E^. For comparison between 
experimental results and theory this distribution function has been 
utilized, with corrections to account for the finite energy width of 
the incident electron beam. This treatment is described in Section IV. 
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III. EXPERIMENTAL PROCEDURE 



The experimental arrangement used was similar to that described 
by Breuer [6*], Miller [ 12 1 , and Goodwin [131. The LINAC of the Naval 
Postgraduate School was used to accelerate electrons to energies from 
54 to 97 MeV. These electrons were elastically scattered at 90° from 
a thin (0.076 mm) aluminum scattering foil, thence passing through the 
absorber, and finally energy-analyzed by a 120° magnetic spectrometer 
described by Kenaston, Luke and Sones [ll]. A figure of this arrange- 
ment is given in Refs. [ 12 1 and [13]. 

According to Miller [12], the results of energy loss experiments 
with aluminum appeared to be quite sensitive to geometrical considera- 
tions, In his work data were collected with the absorbing material 
positioned approximately 22 cm from the scattering foil. The results 
with this arrangement showed a general lack of agreement with theoreti- 
cally predicted half-widths. Only when the absorber was placed approxi- 
mately 3 cm from the scattering foil did the experimental half-widths 
show reasonable agreement with theoretical prediction and the results 
of Breuer [6], 

Consequently, an experiment was conducted to determine the optimum 
positioning of the absorber with respect to the scattering foil. A 
device was fabricated whereby the absorber could be positioned at 
various distances from the scattering foil. For details of the experi- 
mental arrangement and a diagram of the apparatus, see Figs. 1 and 2. 

A portion of Miller's experiment was repeated, using an incident energy 

2 

beam of 53.8 MeV and an aluminum absorber of thickness 1.441 gm/cm . 

The separation was varied from 3 to 20 cm and the energy distribution 



14 



measured. The results of these measurements as well as the theoretical 



prediction are shown in Fig. 3. It was determined that the separation 
had no effect on the location of the peak of the energy loss curve, 
but that the half-width increased as the distance increased. Agreement 
with theory as well as previous experimental data was attained for a 
radial separation < 6 cm. Thus in the remainder of this work the ab- 
sorbing layers were positioned approximately 3 cm from the scattering 
foil. 

Because the incident electron beam scattered into the absorbing 
material is not monoenergeti c , the energy distribution of the electron 
beam was measured both before and after passing through the absorbing 
materials. Various thicknesses of each material were positioned into 
the beam by a remotely-controlled device. This enabled changing of 
absorber thickness without turning off the accelerator and thus poss- 
ibly altering the character of the electron beam. The data represent 
the number of electrons detected by a coincidence counting system at 
the exit of the magnetic spectrometer. A downstream Secondary Emis- 
sion Monitor was used ; as a standard for normalization purposes, in 
that each data point corresponds to a given predetermined integrated 
current. This corresponded to a certain number of electrons passing 
through the absorbing material. 
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IV. TREATMENT OF DATA 



It was first necessary to correct all data taken for background, 
counting loss and accidentals. Background was measured by taking 
counts with the scattering foil removed from the path of the incoming 
beam and integrating to the appropriate current. Counting loss 
corrections were applied using the formula: 

n = n(l + Kn) (11) 

corr 

where K is a characteristic constant of the counting system - in this 
case .003 - and n is the count rate. Accidental coincidence counts were 
measured by means of a 63 nsec delay line, and appropriately deducted. 
Because of low background and the counteracting characteristic of 
counting loss and accidental coincidences these corrections had neg- 
ligible effect on the data. 

Since the theory of Blunck and Westphal is based on a monoenergetic 
incident beam, the energy spread of the incident electrons must be 
taken into consideration. Such an energy spread will increase the 
half-width, and because of the asymmetry of the energy distribution 
will increase the most probable energy loss. To account for these 
effects the IBM 360/67 computer of the Naval Postgraduate School was 
used to unfold the measured incident energy distribution into the 
theory. The method is described in detail by Miller [121 and Goodwin 
[13]. For a computer program to accomplish this unfolding, see Appen- 
dix C. 

In the comparison of the data with theory, the measured incident 
energy distribution was unfolded as described above, properly normal- 
ized to the experimental data and plotted. From these plots the 
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theoretical half-widths and most probable energy losses were extracted 
The experimental data were superimposed on the theoretical curves and 
the measurable parameters compared. The results for various incident 
energies, absorber materials, and thicknesses are shown in Fig. 4 
through 23. The solid lines represent the theoretical energy distri- 
butions, while the crosses indicate experimental data points. To 
avoid clutter on the figures approximately 257> of the data points have 
been omitted. A representative sampling of error bars is included. 
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V. RESULTS AND OBSERVATIONS 



The theoretically predicted and experimentally measured values of 
half-widths and most probable energy losses for the various incident 
energies, absorber materials, and thicknesses are shown in Tables I, 

II, and III. In addition to the experimental data and the Blunck and 
Westphal predictions, the theoretical value of the most probable energy 
loss due to ionization, S, is listed for comparative purposes. The 
data for an incident energy of approximately 75 MeV for all three 

materials of various thicknesses are plotted in Figs. 24 and 25. 

2 

For the thickest lead absorber (t = 2.825 gm/cm ) at 91.37 MeV and 

2 

the copper (t = 5.726 gm/cm ) at 53.63 MeV a half-width was not ob- 
tained because of the shape of the low energy tail. In an attempt to 
get a half-width the energy spectrum of the electrons was examined 
down to about 5 MeV, the lower limit of reliability of the magnetic 
spectrometer. It was found that for thick, absorbers secondary effects 
(probably due to shower production) cause a secondary peak of low 
energy electrons. This rise began in the vicinity of 30 MeV and peaked 
at about 15 MeV, approximately independent of incident energy and ab- 
sorber material. It was noted that the amplitude increased with target 
thickness, and in the case of the thickest absorbers was as much as 
ten times greater than the primary peak. Measurement of the presence 
of positrons was made by reversing the polarity of the magnetic 
spectrometer. The energy spectrum of the positrons at low energy 
followed that of the electrons, but the intensity was much lower. This 
indicated that the lower energy tail was due in part to pair production 
resulting from the interaction of brems trahlung with the absorber 
material, inasmuch as only a negligible number of positrons were 
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detected when the absorber was removed. Another possible contri- 
bution of electrons to the low energy tail is electron-electron 
scattering . 

From Table I it may be seen that for lead the experimental results 

show agreement with theory in the most probable energy loss for all 

2 

energies and thicknesses up to 2.825 gm/cm . There is agreement in 

2 

the half-width for thicknesses up to 2.118 gms/cm , but at greater 

thickness the experimental half-width is smaller than that predicted 

by theory. This is in contrast with previous works with materials 

of lower Z where for thick absorbers the experimental half-widths 

are systematically greater than theory predicts. 

2 

Table II indicates that for a 5.574 gm/cm aluminum absorber the 
most probable energy loss data agrees with theory within the error 
limits of the experiment. However, at an incident energy of both 
53.44 and 74.47 MeV the experimental half-width is significantly 
greater than theory predicts. 

For copper, again agreement between experimental data and theory 

2 

for the most probable energy loss was observed up to 5.726 gm/cm . 

There appears to be agreement in the half-width for a thickness of 

2 2 
4.295 gm/cm , but the results of the experiments using the 5.726 gm/cm 

absorber are inconclusive owing to the large degree of uncertainty. 

It would appear in general that the Blunck and Westphal theory for 
the energy loss distribution of electrons is valid for predicting the 
most probable energy loss for aluminum, copper, and lead for incident 
energies in the range 54 to 97MeV. The theoretical prediction of half- 
width for thick absorbers seems best for copper, but perhaps is sys- 
tematically too large for materials of higher Z and too small for 
those of lower Z. 
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In light of the foregoing further work is justified with materials 
of Low and intermediate atomic number. Furthermore, a closer exami- 
nation of the characteristics of the low energy tail and secondary 
peak is warranted. 
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APPENDIX A 



TABLES 



TABLE I. Energy Loss Distribution Characteristics of Lead 





Q (MeV) 
P 


HW (MeV) 


E. (MeV) 


2 

t (gm/ cm ) 


B&W 


S 


Experiment 


B6W 

Folded 


Experiment 


53.85 


0.706 


0.83 


0.69 


0.81 + .05 


0.61 


0. 60 + . 06 




1.412 


1.41 


1.43 


1.45 + .08 


1.21 


1.23 + .12 




2.118 


2.65 


2.20 


2.64 + .10 


2.51 


2.65 + .25 




2.825 


3.75 


2.98 


3.73 + .15 


7.44 


3.00 + .30 


74.74 


0.706 


0.85 


0.69 


0.82 + .05 


0.83 


0.82 + .05 




1.412 


1.80 


1.45 


1.80 + .10 


1.40 


1.47 + .07 




2.118 


2.70 


2.22 


2.59 + .12 


2.50 


2.24 + .19 




2.825 


3.90 


3.01 


3.80 + .20 


8.78 


6.03 + 1.0 


91.37 


0.706 


0.92 


0.70 


0.87 + .07 


0.89 


0.78 + .08 




1.412 


1.84 


1.45 


1.77 + .10 


1.54 


1.52 + .19 




2.118 


2.84 


2.23 


2.75 + .15 


2.85 


2.74 + .64 




2.825 


4.04 


3.02 


3.88 + .20 


* 


* 



* No half-width obtained. 
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TABLE II. Energy Loss Distribution Characteristics of Aluminum 





Q (MeV) 
P 


HW (MeV) 


E. (MeV) 
1 


tCgm/cm 2 ) 


B&W 


S 


Experiment 


B&W 

Folded 


Experiment 


53.44 


5.574 


9.22 


8.41 


9.24 + .13 


3.65 


4.60 + .30 


74.47 


5.574 


9.13 


8.43 


9.33 + .25 


4.00 


5.7 + 1.2 



TABLE III. Energy Loss Distribution Characteristics of Copper 





Q (MeV) 
P 


HW (MeV ) 


E (MeV) 


2 

t (gin/ cm ' 


B&W 


S 


Experiment 


B6W 

Folded 


Experiment 


53.63 


4.295 


6.28 


5.56 


6.23 + .13 


4.30 


4.25 + 1.20 




5.726 


8.90 


7.52 


8.83 + .20 


* 


* 


76.02 


4.295 


6.29 


5.57 


6.28 + .13 


4.59 


4.55 + .50 




5.726 


9.13 


7.54 


8.38 + .25 


9.20 


10.1 + 2.0 


94.07 


4.295 


6.30 


5.58 


6.40 + .20 


4.61 


~ 4 


96.66 


5.726 


9.30 


7.55 


9.2 + .5 


15.55 


16.4 + 4.2 



*No half-width obtained. 
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APPENDIX B 




FIGURE 1. Experimental Arrangement for Varying Foil-to-Absorber 
Distance in Beam Line. 
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FIGURE 2. Device to Vary Radial Distance 
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FIGURE 3. Half-Width of 1.441 gm/cm A1 Absorber vs Radial Distance 
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FIGURE 4. Pb, E. = 53.85 MeV, t = 0.706 gm/cm 2 
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FIGURE 6. Pb, E. = 53.85 MeV, t = 2.118 gra/cm 2 
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FIGURE 7. Pb, E =53.85 MeV, t = 2.825 gm/cm 
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FIGURE 9. Pb, E = 74.74 MeV, t = 1.412 gm/cm 2 
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FIGURE 12. Pb, E. = 91.37 MeV, t = 0.706 gm/cm 
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FIGURE 13. Pb, E = 91.37 MeV, t = 1.412 gm/cm 2 
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FIGURE 14. Pb, E = 91.37 MeV, t = 2.118 gm/ctn 2 
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FIGURE L6 . Al, E. = 53.44 MeV, 



= 5.574 gm/cm 
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FIGURE 17. Al, E. = 74.47 MeV, t = 5.574 gm/cm 2 
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FIGURE 20. Cu, E i = 76.02 MeV, t = 4.295 gm/cm 2 
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FIGURE 21. Cu, E. = 76.02 MeV, t = 5.726 gm/cm 2 
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FIGURE 23. Cu, = 96.66 MeV, t = 5.726 gm/cm 2 
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FIGURE 24. Most Probable Energy Loss at E^ « 75 MeV 
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APPENDIX C 



dpO^rah Fro COMPUTING ENERGY LOSS D 1ST R I B'JT I CM 



IMPLICIT REAL*8 (A-H,L,Q-Z) 

DIMENSION E(IOO). C7(100)» EXI200), CT(2CO) 
CDMMON/tVL/EI ,Q,QB ,AR ,7 ,ALR,BETSO 
EXTERNAL 
W°I T E (6.B8) 

po F opm A t ( * 1 ’ ,11X. ’El ’ ,7X, »BETA**2’ ,13X, • AF • , l IX, • ALPHA 
I R ’ ,°X, ’OBAR • » 1 2 X , ’DPI * »////) 

READ N- NUMBER OF I NFUT ENERGY BINS £ M-NUMBER OF OUTPUT 
ENERGY VALUES 

sPAD ( 5,10) N , M 

10 PORM^T (2110) 

R E A n ( c ,ll) R 

11 FORMAT (F10.5) 

READ Z-ATOMIC NR, PHO-DENSTTY, A-ATOMIC WT , £ THE PARA- 
METERS OF THE MATERIAL FROM PAGE 514, STFRNHEl MER PAPER 

RPAn (5.12) 7 ,PHO, A, AS, B» CS , SMALA, EM, XI 

12 FORMAT (6F10.4,/,3F10.4) 

READ 7 F° 0 FEAK ENERGIES AND COUNTS FROM HIST OGR A M 

READ (5,13) (F( J) ,C7(J) ,J=i,N) 

13 phrmaT ( F 1 G. 2 ,F 1 0. 0) 

D n ICC 1=1, M 

ICC CT( I ) = C.C 

F X ( 1 ) IS THE HIGHEST OUTPUT ENERGY VALUE DESIRED 
EX( 1 ) = SC.0C 

f'LTE IS THF SEPARATION RE T WFEN ENERGY OUTPUT VALUES 



DL T E = • I 

rOMDUTE ALPHA *R 

ALR = . CC 14 * 7 #* 2 *RH 0 *R/A*( 4 . / 3 . *DLOG ( 1 8 3 . / ?** ( 1 .73 . ) ) 

l + l./<3.) 

T = P HO*P 
DO 2 C 0 J*l,N 
E I=E ( J) 

COMPIJ T E B F T A * *2 

BE T SO=l.-. 511 ** 2 /(El +. 511) **2 
CQMPU T E A *R 

AR=. 154*7*RHC*P /( A*BETSQ) 

COMPUTE P/(M*C) 

X=DS 0 RT(EI** 2 + 1 . 022 *EI ) / • 5 1 1 
COMPlJTf q-bap 

OB=AS*T/BETSO*( B+. 43 +DL 0 G( El ) -BETSO-CS -S MAL A* (XI 
1 -DLOG 1 C( X) )**FM) 

COMP U T F Q PROBABLE I 0 NI 7 .ATI 0 N 

OPI=A S*T /BE TSO*(B-H. 06 + DL 0 G( A S*T/BETSO) -BETSO-CS 
1 - SMA LA * ( Xl-DLOGIO(X) ) ♦*E M ) 

DO 300 1 = 1 , M 
Q=F I-E X( I ) 
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rx( i+n=Fxn i-clte 
ip*o.l=.c.O) or T c 3CG 

lirv \ DC TH 17 I\ T FDM£PT A T F LI M I ~8 C f JNJTcpq at j pm 

U=. CCC0CC1 *0 

V=. 2 c *o 

cm 0003210 . o,u,w,fi ) 

CALL 0C0?2H!,V,W,F2) 

CUL 000321 V,Q,W,F3) 

F = F l + R 2 + F 3 

C T { I >=C T < I H-F*f7( J) 

IOC C^NTJMUF 

W^ITF (6,23) ET ,RFT$0,A R,AL3,0* ,OPI 

20 Ff)PM>T (C 15 . 2 ,F 1 *. 9 , 4 F 15.4 ,/) 

2CC CONTJNUF 

WO T T ? ( 6 , 7 7 ) 

7-7 chpm»T (• 1 1 ,6X , * c NEPCY • ,3X , 'CCUNTS *,//// ) 
WRTtf (A,21) (FX(I),CT(I), I=1,M) 

21 format (F11.2,F15.0) 

RTpp 

END 



FUNCTION TO C OM D UTE THE ENERGY LOSS OWE TO ICNIZATIPN 



PEAL FUNCTION WI *8(EI ,Q , GB , A R , 7 , RFTS Q) 
implicit R=AL*3 IA-H,L,0-Z) 

CTORE THE RLUNC* r. LE I C FG8 NG CONST ANTS 

REAL* 0 S«4) /C. 174,. CE 8,. 01°,. 00 7/, G4MM$ (/,)/l .3,2.0, 3.0 
l, 3.0/, LAMBDA (4) /O. 0,3. 0,6. 5, 11.0/ 

STOPF THE IONIZATION POTENTIALS AND ELECTRON LEVEL* 7 OF 
THE MATERIAL 

pfAL*8 I m ( 6) /. 0878 9, .01432, .002 01, .00038 , .OCOO^o, 

1 .00003? / , NM ( 6) /2. ,8. ,18. ,32. ,13. ,4. / 

W IPN = G * 0 
RR=C.C 

CALCULATE la mroA 

ALA M = ( O-OB ) / A R+D LOG (El /AR ) -1 .116 
CALCULATE B **2 
DO ICC 1*1,6 

R=!M( I )*NV(I )*DL0G<2. *FI /( TM(I)*(1.-BE T S0) ) ) 

ICO SB = SB +B 

BS2=3.*$B/( 7*AP) 
no 400 K = 1 ,4 

Fp = ( ( A LA ^-LAMBDA IK) ) ««2) /( GAMMS (K) **2 + BR 2 ) 

XW=0F XP ( — c ° ) 

TF( XW.LF. 1. C0-60) XW=0.0 

WIN*( S( K ) *(C amm.S( K) / AR) /OSORT (GAMM* 7 ( k ) **? + RR 2 )) * XW 
a 00 WIONNMCN+WI N 
W I = W 1 0 M 
DETypM 
END 



FUNC T I CN TO COVpUTF THE E NE°GY LOSS DUE TO RADI AT ION 



R FA L FUNCTION VS *8 < F I , X . AL R ) 

IMPLICIT REU *8 ( A-H , L , 0-7 ) 

W* = ALR/DGA «MA ( ALP+ 1 . ) *(X/FT )**ALF/X 
RETURN 

FN|p 
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FIJNCTTCN TQ COMPUTE TOTAL C NFRGY LOSS 



REAL FUNCTION W*R(X) 

IMPLICIT REAL*« (A-H , LtO-Z ) 

Common / tVL /El ,0,QB ,AP ,7 ,ALR,8ETSQ 
W=WT( El t Q-X,OB,AP,Z,BET$Q) *WS(EI ,X ,ALR) 
RE T UPN 
ENO 



SUBROUTINE DQG32(XL»XU,FCT,Y) 



OUROCFE 

TC COMPUTE INTFGRAL(FCTIX) 

USAGF 

CALL OCG32 ( XL , XU, FCT.Y > 

PARAMETER FC T PEOUIRFS AN EXTERNAL STATEMENT 

DESCRIPTION OF PARAMETERS 

XL - DOUBLE PRECISION LOWER BOUND 

XU - DOUBLE PRECISION UPPER BOUND 

RCT - FUNCTION TO BE INTEGRATED 
Y - THE RESULTING DOUBLE PRECISION VALUE 

ME THOD 

EVALUATION IS DONE BY 32-RCINT GAUSS QUADRATURE 

DOUBLE PRECISION X L , XU , Y , A , B , C , FCT 
A=. 5DO*( XU+XL) 

B = XU- XL 

C=. 49 8631930924 74078D0*B 

Y=. 350930 5004735C4 830-2 *( FCT ( A+C ) +FCT ( A-C ) ) 

C = . 49 28C 5 755772 6341 7D0*B 

Y = Y+. 81 37 19736 5*- 52 83 5 0-2 *( FCT(A+C) + FCT ( A-C ) ) 

C=. 4 823P1 127793753 22 D0*8 

Y=Y+. 1265 603265463 10300-1 *( FCT (A+C) +FCT (A-C) ) 

C = . 4674530379688698400 *8 

Y=Y+. 171369314* 65 1071 ?D-1*( FCT (A+C) +FCT (A-C) ) 

C =. 44 fl 1605778830260600*8 

Y=Y+. 2141794901 111 33400-1*1 FCT I A+C ) +FCT (A-C) ) 

C=. 424 68 380 6 6662849900*8 

Y=Y+. 25* 9RC2963 11 880880-1*1 FCT I A+C)+FCT( A-C) ) 

C=. 397241 897583 071 20D0*B 

Y = Y+. 25 34204673 92677740-1*1 FCT I A+C) +FCT (A-C) ) 

C=. 366 C° 10 59370144 9400 *8 

Y = Y+. 329111113881 8092 3 D-l * ( FCT ( A+C ) +FCT ( A-C ) ) 

C=. 33 15221 334651 C760D0*B 

Y= Y+. 36 17 28970544242530-1 *( FCT ( A+C ) +<=CT ( A-C ) ) 

C = . 2° 3 857878620381 1600*B 

Y=Y+. 3 C C969478935351 530-1* (FCTIA+C) +FCT ( A-C ) ) 

C = . 2534499544661147000*8 

Y=Y+. 416559621134733 78 D-1*(FCT( A+C) +FCT ( A-C ) ) 

C. = . 210 675 6380653176700*8 

Y=Y + . 4 38260465 02201 9060-1 *( FCT ( A+CJ+FCT (A-C) ) 

C=. 1659343011410638200*8 

Y=Y+.455P69393*7881942D-1*(FCT(A+C)+FCT (A-C) ) 

C=. 119643681 126 068540 0*B 

Y= Y+ . 4 6922 19954040228 3 D-l * ( FCT ( A+C ) +FC T ( A-C ) ) 

C=. 722 3 5980791 3982 50- 1*8 

Y=Y + .4 781936C03 563 743 00-1*(FCT( A+C)+FCT (A-C) ) 

C = . 24 1 538 328438691 580-1 *B 

Y=B ♦ ( Y + . 4827C044257363 9000-1 *(FCT(A+C)+FCT (A-C) ) ) 

RETURN 

FNP 
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